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We analyze the requirements to test some of the most paradigmatic collapse models with a protocol 
that prepares quantum superpositions of massive objects. This consists of coherently expanding 
the wave function of a ground-state-cooled mechanical resonator, performing a squared position 
measurement that acts as a double slit, and observing interference after further evolution. The 
analysis is performed in a general framework and takes into account only unavoidable sources of 
decoherence: blackbody radiation and scattering of environmental particles. We also discuss the 
limitations imposed by the experimental implementation of this protocol using cavity quantum 
optomechanics with levitating dielectric nanospheres. 



I. INTRODUCTION 

In the last decades, seminal experiments have demon- 
strated that massive objects can be prepared in spa- 
tial superpositions of the order of its size. This has 
been realized with electrons [1], neutrons [2], atoms and 
dimers [3] , small van der Waals clusters [4] , fullerenes [5] , 
and even with organic molecules containing up to 400 
atoms [6]. These experiments are designed to observe 
the interference of matter- waves after passing, in essence, 
through a Young's double slit. The possibility of observ- 
ing these quantum pheonomena with yet larger objects is 
extremely challenging. This is due to the great quantum 
control and isolation from the environment that these 
experiments require. 

More recently, the field of cavity quantum electro- 
and optomechanics [7-11] has opened the pathway to 
bring much more massive objects to the quantum regime, 
namely, objects containing billions of atoms, thereby im- 
proving the previous benchmark by many orders of mag- 
nitude. This allows to explore the physics of a completely 
new parameter regime. A first step towards this direc- 
tion has been realized in Refs. [12-14], where ground 
state cooling of mechanical resonators at the nano- and 
microscale has been achieved. Additionally, various re- 
searchers have proposed to exploit the coherent cou- 
pling of the mechanical resonator with single photons 
or qubits to create quantum superpositions, see for in- 
stance [15, 16]. In these proposals, the superposition 
of the mechanical motion state is, typically, of the form 
|0) + |1), where |0) and |1) are respectively the ground 
state and the first excited state of the harmonic poten- 
tial. In these states, the position is dclocalized over dis- 
tances of the order of the zero point motion, i.e. xq = 
y?i/(2mu), where m is the mass of the object and ui the 
frequency of the harmonic potential. Within the mega- 
hertz regime, objects containing n at atoms are delocal- 
ized over distances of the order of lO -7 ??.^ 1 ^ 2 m, which is 
subatomic for objects containing billions of atoms. This 
is in contrast with matter- wave experiments, where, de- 
spite the fact that objects have "only" hundreds of atoms, 
they can be delocalizcd over distances larger than their 
size. 

Remarkably, these experiments might be applied to 



the service of a very fundamental goal, namely, the ex- 
ploration of the limits of quantum mechanics predicted 
by several collapse models [17-26]. The common idea of 
these models is the conjecture that the Schrodingcr equa- 
tion is an approximation of a more fundamental equation, 
which breaks down when objects above a critical mass 
are delocalized over a critical distance. This prediction 
is very difficult to confront because of the following ar- 
gument: standard decoherence [27, 28], described within 
quantum mechanics, also predicts the impossibility to de- 
localize large objects due to the interaction with the en- 
vironment; thus, this masks the effects of collapse mod- 
els. This poses a major challenge to corroborate collapse 
models, as the effects predicted by these must stand alone 
from decoherence processes and be exposed to potential 
falsification. This leads to the central questions of this 
work: How challenging is it to test collapse models while 
also taking into account unavoidable sources of decoher- 
ence? Is it preferable to have small objects delocalized 
over large distances, as in matter-wave experiments, or 
rather large objects delocalized over small distances, as 
in experiments with mechanical resonators? 

The aim of this paper is to address the latter 
questions by analyzing a prototypical experiment that 
bridges approaches from quantum-mechanical-resonators 
and matter-wave interferometry. This experiment re- 
lies, on the one hand, on techniques of cavity electro- 
optomcchanics to prepare a mechanical resonator in the 
ground state of its harmonic potential. On the other 
hand, the experiment mimics matter-wave interferome- 
try, as the ground-state-cooled mechanical resonator is 
released from the harmonic trap in such a way that it 
coherently delocalizes over distances much larger than 
its zero point motion xq. A subsequent measurement of 
the squared position, which is to be realized using tech- 
niques of quantum-mechanical resonators, collapses the 
state into a superposition of different spatial locations, 
thereby acting as a Young's double slit. Finally, the sub- 
sequent free evolution generates an interference pattern. 
We remark that the implementation of this experiment 
using cavity quantum optomechanics with optically lev- 
itating dielectric nanospheres [16, 29-31] has been re- 
cently proposed in [32] . The present article analyses this 
proposal with a broader scope, namely, it studies the ef- 
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feet of some of the most paradigmatic collapse models 
together with unavoidable sources of decoherence. This 
allows us to obtain the environmental conditions, masses 
of the objects, and derealization distances where col- 
lapse models can be falsified. These conditions are gen- 
eral and will be common to any physical implementation 
of the protocol. 

This article is organized as follows: in Sec. II we in- 
troduce and analyze the quantum-mechanical resonator 
double slit experiment in a general fashion, neglecting 
decoherence and without specifying the experimental im- 
plementation. The effect of unavoidable sources of deco- 
herence, such as blackbody radiation and scattering of 
environmental particles, will be the subject of Sec. III. 
The effects of several collapse models in this experiment 
are discussed in Sec. IV, where we also obtain the param- 
eter regime needed to confront them. Finally, in Sec. V, 
we study the restrictions imposed by an implementation 
of this experiment using cavity optomcchanics with op- 
tically levitating dielectric nanospheres. We draw our 
conclusions and provide further directions in Sec. VI. 



II. MECHANICAL RESONATOR 
INTERFERENCE IN A DOUBLE SLIT 

In this section we analyze a protocol that merges tech- 
niques and insights from quantum-mechanical resonators 
and matter-wave intcrferomctry; we call it MEchanical 
Resonator Interference in a Double slit (MERID). We 
analyze it without taking into account standard decoher- 
ence (c/. Sec. Ill) and without specifying its experimental 
implementation (c/. Sec. V). The MERID is realized by 
applying the following steps, see Fig. 1: 

1. Prepare a mechanical resonator. For instance, trap 
a massive object of mass to, which is typically a 
sphere, into an harmonic potential with trap fre- 
quency UJ. 

2. Cool the center-of-mass along one direction, say x, 
to, ideally, the ground state of the harmonic poten- 
tial. 

3. Switch off the harmonic trap and let the wave func- 
tion expand freely during some "time of flight" t\. 

4. Perform a measurement of x 2 , that is, of the 
squared position of the cooled degree of freedom. 
This measurement acts as a Young's double slit 
since, given the outcome x 2 , the state collapses into 
a superposition of being at +x and at —x. The 
mechanical resonator is thus prepared in a spatial 
superposition separated by a distance d = 2\x\. 

5. Let the state evolve freely during a second time of 
flight t 2 . 

6. Perform a measurement of the center-of-mass posi- 
tion x. 




FIG. 1: (color online) Schematic illustration of MERID. a) A 
sphere of mass m is harmonically trapped, with frequency w, 
and cooled into the ground state. The zero point motion is 
given by xq. b) The trap is switched off and the wave function 
expands freely during some time of flight ti. At this time, the 
width of the wave function is given by a. c) A squared posi- 
tion measurement is performed such that the wave function 
collapses into a superposition of two wave packets of size aa 
separated by a distance d. Both aa and d depend on the mea- 
surement outcome, d) The superposition state evolves freely 
during a second time of flight £2- An interference pattern is 
formed with peaks separated by xj. 



Repeat the experiment and collect the data for each 
double slit distance d, corresponding to the result 
of the squared position measurement. An interfer- 
ence pattern in the final position measurement is 
unveiled for each d. 



A. Steps 1 and 2: Cooled initial state 

These steps consist in preparing the object's the center- 
of-mass motion along x in the ground state of an har- 
monic potential with a trapping frequency lo, see Fig. la. 
In the ideal case, the wave function is given by 
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where xq = \f h/ (2tocj) is the ground state size and m is 
the mass of the object. In realistic situations, the initial 
state is given by a thermal state with mean occupation 
number n = (cxp[/3fto;] — 1) _1 (where P~ l = fc&T, kb is the 
Boltzmann constant, and T the effective one dimensional 
center-of-mass temperature) , which can be written in the 
Fock basis as 
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This state has the following moments (x (0)) = (2n + 
l)x 2 , (f(0)) = (2n+l)H 2 /(4x 2 ), and ([i(0),p(0)] + ) = 0. 
We do not discuss here how cooling is performed experi- 
mentally, see however [16, 29-31, 33] for optomechanical 
cooling techniques [34-36] applied to optically levitating 
nanosphcres. 



B. Step 3: Expansion 

This step (see Fig. lb) consists in switching off the 
trap and letting the wave function evolve freely, that 
is, it evolves with the unitary time evolution Uo(t) = 
exp[— iH t/ti\, where H — p 2 /(2m). Considering the 
initial state to be the pure ground state, the state after 
some time t\ is given by 
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where a 2 = x\{\ + t\uj 2 ) is the size of the expanded 
wavcfunction, and <^to£ = uti/4 is the global phase accu- 
mulated during the free evolution. 




FIG. 2: (color online) IV^)! 2 = l(a#)| 2 > see Eq. (4), is plot- 
ted for d = a/2 and measurement strength x = 6 (dotted 
gray), x = 10 (dashed red), and x — 20 (solid blue). 



characterizes the strength of the measurement, and re- 
lates d with ad as 



C. Step 4: Double slit 

In this step a squared position measurement of the 
state at time t = t\ is performed, see Eq. (3), such that 
the state collapses into 
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The measurement operator, Aid, is assumed to have the 
following form 

M d =exp [\cj) ds (x/a) 2 ] x 
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This measurement has the potential to prepare a quan- 
tum superposition of Gaussian wavcfunctions of width 
o d separated by a distance d, with an added global phase 
that we discuss below. The state presents a well- 
resolved spatial superposition provided that d > 2a d- 
Also, one requires ^/8a > d in order to have a non- 
negligible probability to obtain the result d, that is, in 
order to ensure that |(d/2|L> (ii)|<i/2)| 2 /|(0|L>o(ti)|0)| 2 > 
exp[— 1]. We have summarized in Table I all the condi- 
tions required to successfully realize MERID that will be 
obtained throughout the article. These two obtained here 
are included as conditions (i) and (ii) by using the defi- 
nition given below in Eq. (6). Motivated by the optome- 
chanical implementation of this measurement (see Sec. V 
where wc derive Eq. (5)), let us define a dimensionlcss pa- 
rameter independent of the measurement result, which 
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For a given outcome d, the larger the value of x, the more 
resolved the superposition. Figure 2 shows the position 
probability distribution of the state of Eq. (4) with d = 
a/2 for different measurement strengths x- 

Finally, note that a global phase </>ds is added during 
the measurement. This phase, as well as the one accumu- 
lated during the time of flight, to f in Eq. (3), plays an 
important role. The condition |(/>d s + 0tofM 2 /(4c 2 ) <C 1 
needs to be fulfilled in order to build the interference of 
the two wave packets centered at x = d/2 and x = —d/2. 
This can be shown by analyzing (p\ip), that is Eq. (4) in 
momentum space, for different global phases, see Fig. 3. 
Recall that within free evolution, the probability momen- 
tum distribution of a wave function at t\, has the same 
form as the position probability distribution at much lat- 
ter times since x{t\ + t^) ~ j>{t\)ti/m. More intuitively, 
the global phase adds some momentum to the wave pack- 
ets. Depending on the sign of this phase, the two wave 
packets either move apart or towards each other. If the 
momentum given is too large, for the former case they 
will separate with a velocity faster than their expansion 
rate, and thus, they will never overlap. For the latter 
case, they overlap only during the "collision" , however, 
at this time, the wave packets have nearly not expanded 
and the fringes in the interference pattern cannot be re- 
solved, see the discussion below. 



D. Steps 5 and 6: Interference 

These steps consist in measuring the position distribu- 
tion of the state obtained after letting the system evolve 
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FIG. 3: (color online) |i/>(p)| 2 = \(pW\ 2 , see Eq. (4), is plotted 
in arbitrary units for d = a/2, measurement strength \ = 50, 
and a = |</>ds + </>tof| equal to a = (dotted gray), a = 50 
(dashed red), a = 100 (dot-dashed blue), and a = 150 (solid 
purple) . 
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FIG. 4: (color online) The correlation function T(x), see 
Eq. (9), is plotted (solid blue line). The short-limit approx- 
imation T(x) — Ax 2 is also plotted for comparison (dashed 
gray line). 



freely during a second time of flight t2] this reads 
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The state \ipf) presents interference peaks separated 
by a distance Xf = 2ith±ij (md) as long as \<fid s + 
0tofM 2 /(4c 2 ) <C 1. The peaks are clearly visible when the 
two wave packets overlap, that is, when d = t 2 h j '(2a d m ) ■ 
Using Eq. (6) and a 2 ~ x^tlu 2 (valid at t\w 1), one 
obtains an upper bound for t\ given by t\ < sjltixl^; 
this corresponds to condition (iii) in Table I. Another 
condition is given by the requirement to resolve the in- 
terference fringes. Assuming a position resolution of Sx, 
one requires Xf > Sx, which provides an upper bound 
for the slit distance d given by d < 2~Khti j (mSx) ; this 
sets condition (iv) in Table I. The MERID is finished in 
step 7 where the protocol is repeated to obtain a different 
interference pattern for each double slit length d. 



on the expansion time t x (and therefore the superposi- 
tion size d) , as well as to the visibility of the interference 
pattern. This form of decoherence includes the standard 
sources of decoherence mentioned above as well as the 
effect of collapse models, which we discuss in Sec. IV. 



A. Position-localization decoherence 

1. Master equation 

The main feature of a position-localization decoher- 
ence is the exponential decay of position correlations, 
i.e. (x\p(t)\x') oc e~ rt (x\p(0)\x'), where T usually de- 
pends on | a; — x'\. This form is common both to the 
decoherence caused by interaction with the environ- 
ment [27, 28] and to the exotic one caused by collapse 
models [17-26]. The qualitative behavior of this source 
of decoherence is very well described by the following 
master equation given in the position basis 



III. DECOHERENCE 

In the previous Section we obtained conditions (i-iv) 
in Table I for a successful realization of MERID. Note 
that, since t\ and ti are unbounded, conditions (i-iv) 
do in principle allow for the preparation of arbitrarily 
large superpositions. This is the stage when one has to 
take into account the effect of decoherence, which is the 
subject of this Section. We will concentrate on unavoid- 
able sources of decoherence, that is, on the decoherence 
caused by the interaction with environmental massive 
particles (c/. Sec. IIIB1), and the effects of blackbody 
radiation (c/. Sec. IIIB2). We start in Sec. Ill A by an- 
alyzing a general form of decoherence called position lo- 
calization. We derive the limitations that this imposes 



(x\m\x>) = -(x\[p,H]\x') - T(x - x')(x\p(t)\x>), (8) 

where we assume one-dimension for simplicity. The de- 
coherence rate function is defined by 



r(x) = 7 1 — exp 



x 
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This function depends on two parameters: the localiza- 
tion strength 7 > 0, which has dimensions of frequency, 
and the localization distance a > 0, which has dimensions 
of length. The value of these parameters depends on the 
particular source of decoherence, such as gas scattering, 
black-body radiation, or the one given by collapse mod- 
els. This simple master equation captures the important 



5 





ti 


d 


*2 


i) Highly probably outcome 


- 
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ii) Superposition peaks resolved 


- 
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iii) Wave packets overlap 


< v / 2t 2X /co 
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iv) Fringes can be resolved 
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< 2nht 2 /(mSx) 


- 


v) Decoh. expansion if d <C 2a 


Optimal: t max 


< < Cmax 
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vi) Decoh. expansion it a 2> 2a 


< 1/7 






vii) Decoh. interference if d <C 2a 




< V3/(At 3 ) 




viii) Decoh. interference if d ^> 2a 






« 1/7 


ix) Optomechanical implementation 


< niin{ v //c/flro,4flo/rB }/a; 
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TABLE I: . Summary of the restrictions on the expansion time ti, the superposition size d, and the evolution time forming the 
interference pattern ti of the MERID experiment. Conditions (i-iv) are discussed in Sec. II and depend on the measurement 
strength. Conditions (v-viii) depend on the position-localization decoherence and are derived in Sec. III. Finally, condition 
(ix) is given by the optomechanical implementation of MERID, which is the subject of Sec. V. Recall the definitions of t max in 
Eq. (19), of Cmax in Eq. (20), and of £ s (t x ) in Eq. (18). 



feature of position localization decoherence, namely, the 
saturation behavior 



T(x) 



Ax 2 , x < 2a, 
x > 2a, 
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where we have defined the localization parameter A = 
7/(4a 2 ). That is, in the short-distance limit, \x — x'\ <C 
2a, the position correlations decay (ignoring the coher- 
ent evolution given by the Hamiltonian) as {x\p{t)\x') oc 
e -A\x-x | t( x \p((y)\x'), such that the decoherence rate de- 



pends quadratically on 
reads 



In this limit, Eq. (8) 
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The decoherence rate saturates in the long-distance limit 
|ar — £c'| ^> 2a. In this regime, the rate is inde- 
pendent of | a; — x'\ and the position correlations de- 
cay as (x\p(t)\x') oc e~ 7t (x\p(0)\x'). For instance, in 
Sec. IIIB1, we will see that this is the limit where 
the wavelength of the particles impinging the object is 
smaller than the separation of a superposition state, such 
that a single scattering event resolves the position of the 
object and provides which-path information. 

In the rest of the article, our strategy will be to approx- 
imate each source of position localization decoherence by 
the simple master equation of the form (9), obtaining 
the localization strength 7 and distance a. Let us there- 
fore analyze the restrictions that the master equation (9) 
impose on MERID. 



2. Solution of the Master equation 

In MERID, different steps assume the free evolution 
given by the Hamiltonian H = p 2 / (2m). Hence, we need 
to obtain the solution of Eq. (8) for the free dynamics 



case. This is given by [17] 
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x T(p,x - x',t)(x + y\p s (t)\x' +y), 



where p s (t) denotes the evolution of the density matrix 
with the Schrodinger equation only, that is, when r(x) = 
0. The function 
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takes into account the effects of decoherence. We will use 
this solution in the following to obtain different restric- 
tions on MERID. 



3. Free evolution: coherence length 

We have seen that the step 4 of MERID implements 
a double slit. It is crucial to ensure that the squared 
position measurement prepares a quantum superposition 
instead of a statistical mixture. For this to happen, the 
coherence length of the state before the measurement has 
to be larger than the slit separation d. The coherence 
length is obtained by analyzing the decay of the posi- 
tion correlation function C(x,t) = (x/2\p(t)\ — x/2) as a 
function of the distance. 

Let us begin by using the solution Eq. (12) to obtain 
the time evolution of the mean values and moments of 
x and p. It is straightforward to observe that the mean 
values are not perturbed by the position localization de- 
coherence, that is, (x(t)) = (x(t)) s = (x(0)) + t(p(0))/m, 
(p(i)) = (p(t)) s = (p(0)). However, decoherence does 
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modify the time evolution of the second order moments 



(p 2 (t)) = (p 2 (t)) s + 2Ah 2 t, 
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(p 2 (t)) s = 

We remark 



Here, (x 2 (t)) s = (x 2 (0)) + (p 2 (0))t 2 /(2m) 
(p 2 (0)), and ([x(t),p(t)} + ) = 2(p 2 (0))t/m. 
that the extra diffusive term ~ t 3 found in Eq. (14) for 
the position fluctuations is a clear signature of a random 
force without damping, which is the case of the position- 
localization decoherence. It is also interesting to note 
that Eqs. (14) depend only on A and therefore could have 
been obtained with the simpler master equation Eq. (11) 
which is, however, only valid in the short-distance limit. 

The position correlation function C(x,t) can be now 
computed using Eq. (12) and Eq. (14). In particular, 
by taking into account that p s (t) is a Gaussian state, one 
can perform the integration over y in Eq. (12) and obtain 



C(x,t) 



x exp 



(p 2 ) s X 2 + (X 2 ) s p 2 - ([X,p] + ) s xp 

2ft 2 



(15) 



Note that we do not explicitly write the time dependence 
of the moments in order to ease the notation. A simpler 
formula of C(x,t) for the short(long)-distance limit can 
be derived by using the corresponding approximation in 
Tip, x, t). This leads to 
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£ s is obtained by evaluating Eq. (17) with the unitary 
evolution given by the Schrodinger equation alone. Recall 
the definition of a 2 (t) = Xq(1 + t 2 uj 2 ) and that for an 
initial thermal state, (i 2 (0)) = (2n + l)x 2 . While £ s (t) 
increases monotonically in time, has a maximum at 
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FIG. 5: (color online) As an example, the coherence distance 
see Eq. (17), is plotted for a sphere of R = 50 nm at 
a bulk temperature of 200 K taking into account the deco- 
herence given by blackbody radiation, see Sec. Ill B 2. Other 
experimental parameters are taken from Table III. 



See Fig. 5 for a particular example. Notice that the maxi- 
mum coherence distance, according to Eq. (16), depends 
crucially on the saturation distance a of the position- 
localization decoherence. 

As mentioned before, the coherence distance imposes 
some conditions on MERID. In particular, the superpo- 
sition size d has to be smaller than the coherence dis- 
tance, namely, one requires C(d, tx)/C(0, t\) ~ 1 in order 
to prepare a coherent superposition instead of a statisti- 
cal mixture. This gives rise to conditions (v) and (vi) in 
Table I depending on the ratio d/(2a). 



4- Visibility of the interference pattern 

The position-localization decoherence can also compro- 
mise the visibility of the interference pattern in the step 
5 and 6 of MERID. Using Eq. (12), one obtains that 
the time evolution of the position distribution P(x,t) = 
(x\p(t)\x) is given by 



P(x,t) 




dpdx'e* J r (p)e P s (x',t) 



dpe s T(p)P s (p,t), 



(21) 



where P s (p,t) = J dxe~ ipx ' / h P s (x, t)/y/2nh. The posi- 
tion distribution without decoherence P s {x,t) oscillates 
with a wavelength given by Xf = 2nht/ (md), which cor- 
responds to the distances between the interference max- 
ima. Thus, P s (p,t) has peaks at pj = ±2Trh/xf = md/t. 
Hence, the reduction of the interference peaks, which we 
use as a figure of merit for the visibility of the interfer- 
ence pattern, is given by V(t) = J- {2irK/xf) = exp [— iO], 
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Note that 9 w Ad 2 /3 in the limit d < 2a and 9 « 7 
in the limit d S> 2a. Therefore, the requirement 9^2 <SC 
1 establishes the conditions (vii) and (viii) in Table I 
depending on the ratio d/(2a). 



B. Unavoidable sources of standard decoherence 

Let us now focus on the unavoidable decoherence given 
by scattering of air molecules and blackbody radiation. 
Decoherence due to environmental scattering is a well 
studied topic triggered by the work of Joos and Zeh [37]. 
For an extensive review on these topics, we refer the 
reader to the textbooks [27, 28]. Here, we review the 
results needed for our analysis. Localization due to en- 
vironmental scattering is described by a master equation 
of the type 



(x|p(i)|x'> = -<x| [p,H 



|x') -F(x-x')p(x,x'), (23) 



where the decoherence function ^(x — x') depends on the 
distance |x — x'| and can be expressed as [28] 
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The derivation assumes an infinitely massive object and 
the fact that the incoming particles are isotropically dis- 
tributed in space. Here, p{q) denotes the number den- 
sity of incoming particles with magnitude of momen- 
tum equal to q, v(q) = q/m a (v(q) — c) is the veloc- 
ity of massive (massless) particles, |n| = |n'| = 1, and 
/(gn, qn') is the elastic scattering amplitude. For fur- 
ther details, see Chapter 3 of [28]. The behavior of the 
function is very different depending on the ratio between 
|x — x'| and the thermal wavelength of the scattering par- 
ticles Ath- In the long wavelength limit, Ath 3* l x — x '|j 
F(x — x') ~ A|x — x'| 2 , whereas in the short wave- 
length limit, Ath |x — x'|, one obtains the saturation 
of .F(x — x') ~ 7. That is, above some critical distance 
each scattering event resolves the separation |x — x'| and 
thereby provides which path information. This allows 
us to relate qualitatively and quantitatively the master 
equation (23) with the simpler one (8) discussed in the 
previous subsections. This connection, which has been 
discussed previously in [38, 39], is given by the following 
relations 



A th /2 and 7 = A 2 h A. 



(25) 



In the analysis of decoherence due to environmen- 
tal scattering one typically employs the long wavelength 



limit since it always provides upper bounds on decoher- 
ence rates, even when one is in the short wavelength limit. 
This was the case, for instance, in the analysis performed 
in the optomechanical double slit proposal in [32]. As 
shown below, the upper bounds for the case of scattering 
of air molecules were too loose, since one is in the satu- 
ration regime. This yielded the requirement of very low 
pressures. The analysis performed in the following takes 
into account the saturation effect and yields much more 
feasible vacuum conditions. 



1. Air molecules 

The thermal wavelength of a typical air molecule, 
which is assumed to be in thermal equilibrium with 
an environment at temperature T, is given by A^"' = 
2irH/ i l /2irm a ki ) T e = 2a a ; r , where m a is its mass. Using 
m a ~ 28.97 amu and T e ~ 4.5 K, one obtains 2a air ~ 0.15 
ran. The localization parameter associated with scatter- 
ing of air molecules in the long wavelength limit is given 
by [28] 



8V2Trm a vPR z 
3y/3h 2 



(26) 



where v is the mean velocity of the air molecules, P the 
pressure of the gas, and R the radius of the sphere. Thus, 
using Eqs. (25), (26), and the expression of A™, one ob- 
tains 



16ttV2^ PR 2 



V3 



(27) 



In the following, we will consider superpositions which 
are, at least, in the nanometer scale. Therefore, we will 
use the short wavelength limit d 3> 2a to account for 
the decoherence effect of air molecules. The effect of this 
decoherence is shown in Fig. 6, where the coherence time 
7air ano - t ne corresponding coherence distance isd^n) 
are plotted as a function of the diameter of the sphere and 
for different pressures. Note that these quantities define 
the conditions (vi) and (viii) in Table I. In particular, 
for sufficiently low pressures, large superpositions of the 
order of the size of the object are permitted. We remark 
again that in [32] the saturation effect was not taken into 
account, and this gave rise to pressures of 10~ 16 Torr for 
spheres of 40 nm, which turns out to be a very loose upper 
bound when taking into account the saturation effect. 
We will come back to this point in the optomechanical 
implementation of MERID in Sec.V. 



2. Blackbody radiation 



The thermal wavelength for massless particles is given 
by AjJ = 7r 2 / 3 hc/(ki,T e ) = 2abb, which at temperatures 
T ~ 4.5 K takes the value of Ah b ~ 1 mm. In this case the 
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FIG. 6: (color online) Coherence time l/7 a ir (upper panel) 
and the corresponding coherence distance £ s (l/7air) (lower 
panel) taking into account the decoherence of air molecules 
as a function of the sphere's diameter and for environmen- 
tal pressures of P = 10" 11 Torr (solid blue), KT 12 Torr 
(dashed red), 10" 13 Torr (dotted orange), 1CT 14 Torr (dot- 
dashed green), and 10~ 15 Torr (largely dashed purple). Other 
experimental parameters are taken from Table III. In the 
lower panel, the thinner dashed gray line corresponds to the 
line £ s (l/7air) = D. 
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FIG. 7: (color online) t max taking into account blackbody ra- 
diation is plotted as a function of the internal temperature of 
the sphere. We used a sphere of radius of 50 nm although the 
dependence on the radius is negligible. Other experimental 
parameters are taken from Table III. 



dielectric constant, which is assumed to be time indepen- 
dent and relatively constant across the relevant black- 
body spectrum, and Ti is the bulk temperature of the 
object, which might differ from T e . 

From the three contributions, the emission localization 
parameter is usually the dominant one since the internal 
temperature is usually larger than the external one, for 
instance, due to laser absorption during the optical ma- 
nipulation of the sphere. In Fig. 7 we plot the optimal 
time t max , see Eq. (19) as a function of the internal tem- 
perature of the object. The dependence of t max on the 
size of the sphere is negligible. Comparing Fig. 7 with the 
upper panel of Fig. 6, one concludes that, for low pres- 
sure, the decoherence due to blackbody radiation will 
be dominant, specially, when the internal temperature is 
different from the external one which is supposed to be 
cryogenic (a few Kelvins). 



long wavelength limit can be employed since the superpo- 
sitions considered will be always smaller than Recall 
that in this limit the relevant quantity is the localization 
parameter. This parameter has three contributions given 
by scattering, emission, and absorption of thermal pho- 
tons, namely, A b b = A b b, sc + A b b,c + A b b,a, which are 
given by 



A 



bb.se — 



x 8C(9)ci? 6 
9tt 



he 



and 



A 



16ir 5 cR 3 lhT t 



bb,e(a) 



189 



he 



Re 



Im 



£bb - 1 
ebb + 2 



£bb - 1 



£bb 



(28) 



(29) 



We refer the reader to [28, 29] for further details. Here, 
C(x) is the zeta Riemmann function, ebb is the average 



3. Limitations on the superposition size 

Let us now summarize the operational parameter 
regime of MERID taking into account standard sources of 
decoherence. We compute the lower and upper bounds 
for d for allowed values of ti and t% according to Ta- 
ble I. We will use two sets of experimental parameters 
for the pressure, the internal temperature, the measure- 
ment strength, and the common ones given in Table III. 
The first set, which is assumed to be feasible, assumes 
an environmental pressure of P — 10~ 12 Torr, internal 
temperature of the object Tj = 200 K, and measurement 
strength \ = 1000, whereas the challenging set assumes 
P = 10~ 16 Torr, T l = 4.5 K, and X = 10 6 . 

Figure 8 shows the upper bounds for t\ given by 
conditions (iii), (v), and (vi) in Table I, considering 
ti = 0.1/7 a i r . For the feasible set, t\ in the few mil- 
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FIG. 8: (color online) Upper bounds for t\ due to blackbody 
radiation and scattering of air molecules as a function of the 
diameter of the sphere for P = 10" 12 Torr, Ti = 200 K, and 
X = 1000 (P = 10~ 16 Torr, T t = 4.5 K, and \ = 10 6 ) in 
the upper (lower) panel. t max (dashed red), see Eq. (19), is 
condition (v) in Table I, 0.05/7air (dotted brown) is condition 
(vi), and y / 2t2x/u (dotdashed blue) is condition (iii). We 
used t2 = 0.1/7 a i r and the others parameters given in Ta- 
ble III. The shadowed region corresponds to ti fulfilling all 
three conditions. 



liseconds are possible for spheres up to a diameter of 250 
run. For the challenging set, much larger timescales of the 
order of seconds are possible for objects in the microm- 
eter regime. In both cases, pressure is the limiting fac- 
tor for the internal temperatures considered. Note that, 
in both cases, these values correspond to very long co- 
herence times comparing to typical quantum-mechanical 
experiments. 

In Fig. 9 we plot the superposition size that fulfills con- 
ditions (ii), (iv), (v), and (vii) in Table I for ti = 0.1/7 a i r 

and ti = min | y / 2i2x7^-'? ^max, 0.05/7 a i r |, which ensures 

the fulfillment of the restrictions imposed into t\ and ti. 
Hereafter, we will call this plot a a d vs. D" diagram. For 
the feasible set, superpositions larger than the size of the 
sphere are in principle possible for objects of the order of 
50 nm with much lower pressures than those used in [32] . 



FIG. 9: (color online) Superposition distance d as a function 
of the diameter D = 27? of the sphere taking into account 
blackbody radiation and scattering of air molecules with P = 
10" 11 Torr, Ti = 200 K, and X = 1000 (P = 10~ 16 Torr, 
Ti — 4.5 K, and x = 10 6 ) in the upper (lower) panel and 
the others parameters given in Table III. Note the different 
scales for the upper and lower panels. According to Table I, 
we plotted condition (ii) d > a/^/x (solid black), condition 
(iv) d < 2nht2 1 '(mSx) (dotdashed blue), condition (v) d < 
£(ii) (dashed red), and condition (vii) d < ^/3/(Ai2) (dotted 
orange). The shadowed region corresponds to d fulfilling all 
four conditions and thin dashed gray line to d — D. We used 

<2 = 0.1/7air and ti = min j ^2t 2 x/^, tmax, 0.05/7ai r }• 



For the challenging set, larger objects and larger superpo- 
sitions are obtained (note the different scale in the plot). 
It is important to remark that in both cases the limitation 
is given by condition (iv), which reads d < lnhtil (m8x) 
and is related to the resolution of the interference fringes. 
This shows that while coherence times are very large, the 
dynamics resulting in the spreading of the wave function 
are very slow for large masses. This hints at possible 
improvements of MERID using more efficiently the long 
coherence times allowed by the unavoidable sources of 
dccohcrcncc considered here. 
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IV. COLLAPSE MODELS 



nucleon collapse rate, and 



In this Section we discuss the possibility to use 
MERID to test some of the most paradigmatic collapse 
models. Remark that any experimental evidence of these 
would imply a breakdown of the theory of quantum me- 
chanics. In the following we consider the unavoidable 
sources of decoherence discussed in the previous Section 
and. one these grounds, we determine the experimen- 
tal parameters required to falsify a given collapse model 
by the observation of the interference pattern. We note 
that the corroboration of the collapse model is more chal- 
lenging than its falsification, since one must discard that 
standard decoherence, in any of its forms, is responsible 
for the disappearance of the interference pattern. 

We shall not review the extensive literature on collapse 
models; instead, we will focus on some of the most dis- 
cussed ones in the literature, and for each of them, we will 
provide a brief summary of their prediction. In particu- 
lar, we shall express them in a common form, namely as a 
master equation describing position-localization decoher- 
ence; this will allow us to apply the results of Sec. Ill A, 
and thereby to discuss the possibility to test them using 
MERID. 



A. Continuous spontaneous localization 

We start with the Continuous Spontaneous Localiza- 
tion (CSL) model [18, 40], which is the best devel- 
oped collapse model at present [19]. This model builds 
upon the previous works of Ghirardi, Rimini, and Weber 
(GRW) [17], Pcarlc [41], Gisin [42], and it bears some 
similarity to the works of Gisin [43] and Diosi [44, 45]. 
The model is constructed by adding a non-linear stochas- 
tic term to the Schrodingcr equation. This term predicts 
a localization whose strength is directly proportional to 
the mass of the object. At the same time, it is con- 
strained by the fact that the equation must reproduce all 
phenomenology of quantum mechanics for small objects. 
This introduces two phcnomcnological constants that are 
bounded by experimental evidence. 

More specifically, within the CSL model, the master 
equation describing the wavefunction of N particles is 
given in the position representation |xi, . . . , Xjv) = \x) 
by [18, 40, 46] 



(x\lp,H]\x')-r CS L(x,x')(x\p(t)\x'), 



(x\p(t)\x') 
where 



2 " m 

*,3 = 1 



(30) 



(31) 



x - Xj ) + $(x', - x ;j - 2$(x l - x ;j] . 

Here, rap is the mass of a nucleon, 7cg L is the single 



<3>(r) = exp 



An 2 



(32) 



is the localization function with acsL being the local- 
ization distance. Note that for the single nucleon case, 
Eq. (31) reads 



CSL 



(x,x') 



7csl 



1 — exp 



,/|2 



4 «CSL 



(33) 



which has the same form as Eq. (9). The parameters 
Tcsl and acsL are the two phenomenological constants 
of the model. Their value is bounded by both experi- 
mental data and "philosophical" reasons; see [47] for a 
recent discussion. The standard values originally pro- 
posed in the GRW model [17] are acsL = 100 nm and 
Tcsl = 10~ 16 Hz. However, the value of 7csl nas been 
recently reconsidered by Adler and is predicted to be 8 
to 10 orders of magnitude larger [20, 48] than the origi- 
nal one of 10~ 16 Hz, a prediction not yet confronted by 
up-to-date experiments [20]. We will however consider 
here the original values for the sake of comparison. 

The decoherence factor of the CSL model, Eq. (31) 
can be obtained for the center-of-mass of a solid sphere 
of mass m, volume V, and homogeneous mass density. 
We use the results given in Ref. [46] where an analy- 
sis of the CSL model for the free propagation of a solid 
mass is discussed. In this case, the decoherence factor 
(c/. Eq. (31)) takes the form 



r CSL (x,x') = -7cSL?2 



III? 



drdr' 



(34) 



[$(r - r') - $(r - r' + x - x') 



In order to approximate Eq. (34) to (9), we can extract 
the localization parameter in the expression given for the 
free evolution of the position fluctuation, which is given 
by [46] 



(^ W )^(i- 2 (t)) s + 7 H 2/ ^ 2 /aCSL)t3 - 

6m O a CSL 

By comparing it with Eq. (14), one obtains 



ffi^ 7^gj^ 



(35) 



(36) 



(37) 



and has the following limits f(x — > 0) = 1, /(l) = 0.62, 
and f(x — > oo) = 6/x 4 . Thus, recalling that A = 
7/ (4a 2 ), one obtains the collapse rate 



m 2 Aa 2 CSL 
where the function f(x) is given by 

/(*) = 3f 



" 2 












X 2 





7CSL = — 2 7csL./'( i ?/ a CSL). 



(38) 
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FIG. 10: (color online). Different values of the coherence time 
l/(AcsLd 2 ) as a function of the sphere's diameter D and the 
superposition distance d. Other physical parameters, such as 
the density of the sphere, are taken from Table III. 



Note that the rate 7csl grows quadratically with the 
number of nucleons for spheres smaller than 2acsL ■ 

To grasp the strength of the exotic decoherence by this 
model, we plot in in Fig. 10 the value of the coherence 
time 1/(Acsl<^ 2 ) of a superposition of size d <C 2ocsl as a 
function of the sphere's diameter D and the superposition 
distance d. Coherence times of the order of milliseconds 
are obtained for objects of 300 nm and superpositions 
of tens of nanometers. Note that these coherence times 
would be strongly reduced by using the enhancement of 
the localization rate Jqsl predicted by Adler [20, 48]. 



B. Quantum gravity 

Ellis, Mohanty, Mavromatos, and Nanopoulos sug- 
gested in [24, 49] that quantum gravity (QG) can induce 
the collapse of the wavefunction of sufficiently massive 
objects. They argue that the collapse is induced by the 
interaction of the massive object with topologically non- 
trivial spacetime configurations (wormholes) which are 
small compared to physical scales but much larger than 
the Planck scale. From the point of view of quantum 
information theory, this decoherence mechanism can be 
understood as the result of the center-of-mass of the ob- 
ject becoming entangled with some degrees of freedom 
belonging to these wormholes which are unaccessible and 
therefore have to be traced out. Indeed, this provides a 
localization effect analogous to the one induced by the 
interaction with the environment. 



In this model, the quantum wormholes are assumed 
to be in a gaussian state in momentum space, and 
have zero mean momentum with a spread given by 
A ~ ctuq/ (hmp) ~ 10~ 3 m, where mo is the mass 
of a nucleon and mp is the Planck mass. For dis- 
tances smaller than 1/A, the off-diagonal elements of the 
density matrix in position basis decay as (x\p(t)\x') cx 
cxp [— Aq G (x — x') 2 t\ (x\p(0)\x'), where the localization 
parameter is given by 



A 

QG 



h 3 



(39) 



Note that the localization distance is very large since it 
is given by oqg = 1/(2A) ~ 10 3 m, hence, the decoher- 
ence does not saturate within the typical wave function 
sprcadings. The collapse rate for the single nucleon is 
thus given by 7q G = 4oq G Aq G . It is remarkable that this 
model [24, 49], which is based on quantum gravity, con- 
verts the CSL model into a parameter-free model. The 
extension of the single nucleon case to the solid sphere 
can be obtained by retrieving the result given in Eq. (38). 
Based on this, the localization parameter of a solid sphere 
predicted by this model is given by 



A 



QG 



h 3 



(40) 



where we have used /(.R/oqg) = 1 since R <C oqc for 
the spheres considered here. 

As previously, we plot in Fig. 1 1 different values of the 
coherence time 1/(Aqgc? 2 ) as a function of D and d. By 
comparing these results with the ones given by the CSL 
model in Fig. 11, we notice that the decoherence effect is 
slightly stronger, but nonetheless very similar. 



C. von Neumann-Newton equation 

In the last 30 years, many authors have investigated 
the possible role of the Newtonian gravity in the collapse 
of the wavefunction. From these, the independent but 
similar works of Diosi and Penrose (DP) are the most 
famous ones [22, 23, 25, 44, 50-52]. The prediction of 
the model can be castcd into the so-called von Neumann- 
Newton equation, which can be expressed as [23, 44, 50] 

(x|p(t)|xO = 1 (x| [p, H}\x>) r DP (x, x') (x|p(t) |x') , 



where 



(41) 



r . m ?Z g (x,x)-£/ g (x',x') + 2;y g (x,x') 
1dp(x,x)= — . (42) 

Here E/ S (x, x') is the Newtonian interaction between two 
mass densities corresponding to two spheres centered at 
position x and x', respectively. This reads 



C/ 5 (x,x') = -G J 



/(r|x)/(r'|x'h 



drdr\ (43) 
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FIG. 11: (color online). Different values of the coherence 
time l/(AQGd 2 ) as a function of the sphere's diameter D and 
the superposition distance d. Other physical parameters are 
taken from Table III. 
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FIG. 12: (color online) Different values of the coherence time 
l/(ADpd 2 ) as a function of the sphere's diameter D and the 
superposition distance d. The dashed gray line marks the 
limit where the coherence time is not longer valid since, for 
d > 2aup, it would then be given by 1/tdp. Other physical 
parameters are taken from Table III. 



where /(r|x) is the mass density at location r for the 
sphere centered at x. For a rigid homogeneous ball, the 
mass density is uniform and equals / = 3M/(47ri? 3 ). In 
this case, the decoherence function of Eq. (42) depends 
on the relative distance |x— x'| and presents the following 
limits [50] 



DP 



(x) 



Gm 2 /(2R 3 h)x 2 , 
6Gm 2 /(5Rh) 1 



x < R, 
x > R. 



(44) 



The quadratic dependence at short distances allows us to 
identify the localization parameter of the model 



A D p = 



Gm 2 
2R 3 h' 



(45) 



as well as the saturation distance at 2arjp — R. 

The strength of this model is much weaker than the 
CSL and the QG, see Fig. 12 where l/(A D p<i 2 ) is plotted 
as a function of D and d. Coherence times of the order of 
milliseconds are only obtained for objects of few microns 
prepared in superpositions within the micrometer scale. 

Let us point out that the decoherence rate given by 
this collapse model can be strongly enhanced by consid- 
ering the mass density at the microscopic level [23]. For 
instance, by modeling the fine structure beyond the con- 
stant average mass as a conglomerate of identical small 
balls of mass mo and radius ro, the parameters of the 
model arc given by [23] 25dp = '"o an d 



I - ) 



^DP- 



(46) 



Thus, since ro is typically chosen much smaller than R, 
the localization parameter is greatly enhanced. This 
is used, for instance, in the Marshall proposal to test 
the Penrose model using the small delocalization of a 
micro-mirror [15, 53]. However, this choice is contro- 
versial since it does not only convert the parameter free 
model into a one parameter (ro) model, but the von- 
Neumann-Newton equation also becomes divergent for 
point-like particles [23, 52, 54]. This yields unphysical 
results such as the non-conservation of energy, specially 
when r <C 100 nm (see [23, 54]). For this reason, while 
it is not clear how to choose the mass distribution, in 
this article we assume the well-behaved case of a solid 
homogeneous density; this comes at the price of mak- 
ing formidable the possibility of falsifying the model, as 
shown below. 



D. Imprecise space-time 

Finally, we also consider the K-modcl, named after 
Karolyhazy, who introduced one of the first collapse mod- 
els already the 1960's [55, 56]. The model builds upon the 
insight that the sharply determined structure of space- 
time is incompatible with quantum mechanics and gen- 
eral relativity: according to quantum mechanics, the po- 
sition and the velocity of an object cannot have deter- 
ministic values simultaneously, while general relativity 
states that the space-time structure is determined by the 
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TABLE II: Summary of the decoherence parameters predic- 
ted by different collapse models. Recall that, for the CSL 
model we take the original value of 7csl ~ 10 -16 Hz. 

positions and velocities of the masses. 

We base our approach on the article of Frenkel [21] 
who provides a very clear review of the K-model and its 
relation to the CSL model. The prediction of the K- 
modcl in the so-called "no-breathing limit" [21] is given 
by the following master equation 



p(t),H -A K [x,[x,p(t)]] 



(47) 



This corresponds to a position-localization decoherence 
with a localization distance 2a k — > oo, a localization rate 
7x — > 0, and a localization parameter 7^/(4(1^) — > A^. 
The localization parameter for a solid sphere of mass m 
is given by [21] 



2000 F 



1500- 



g iooo- 

500 




500 1000 1500 

D(nm) 



2000 



FIG. 13: (color online) Different values of the coherence time 
l/(AKd 2 ) as a function of the sphere's diameter D and the 
superposition distance d. Other physical parameters are taken 
from Table III. 



A A - = 



8maf ' 



(48) 



where 



«) 3/2 / 



R > ax, 
R < ajc- 



(49) 



Here lp = \J GH/c 3 is the Planck length and lc = H/ (mc) 
the Compton wavelength. Note that this model is also 
parameter- free. 

As shown in Fig. 13, the strength of this model is 
also weaker than the CSL and the QG, and only slightly 
stronger than the DP. Coherence times of the order of 
milliseconds are obtained for spheres with diameter be- 
tween one and two micrometers prepared in superposi- 
tions smaller than one micrometer. 



E. Experimental test 
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FIG. 14: (color online). Localization parameter A for the 
different collapse models (CSL, QG, K, and DP) in units of 
the localization parameter provided by blackbody radiation, 
assuming a bulk temperature of Tj = 4.5 K, as a function of 
the diameter of the sphere. Other parameters are taken from 
Table III. 



Let us now address the possibility to test these col- 
lapse models using MERID. The localization parameter 
and the localization distance of each model are summa- 
rized in Table II. For the sake of comparison, we plot 
the ratio between their localization parameter with the 
localization parameter provided by blackbody radiation 
in Fig. 14, where we have assumed a bulk temperature of 
Ti = 4.5 K (recall Sec. IIIB2). While all collapse models 
provide a stronger localization rate than the blackbody 



radiation, the CSL and QG model are many orders of 
magnitude stronger than the DP and the K-modcl. Ac- 
tually, the standard decoherence given by blackbody ra- 
diation is comparable to the one predicted by the DP and 
K-model for a bulk temperature of 20 K. Since standard 
decoherence will limit the superpositions to be smaller 
than the localization distance, the localization parame- 
ter A is the only relevant parameter. 
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FIG. 15: (color online). Operational parameter regime of 
MERID where the CSL and the QG collapse model can be 
tested. The upper two panels show the d vs. D diagram as in 
Fig. 9 with P = 10" 14 Torr, Tj = 100 K, and X = 1000. The 
green region is the non-overlapping allowed region taking into 
account only standard decoherence and the CSL (first panel) 
or QG (second panel) collapse model. The dot is at d = 30 
nm and D — 100 nm. For this value, the third lower panel 
plots the simulated interference pattern taking into account 
standard decoherence only (dotted gray line), the CSL model 
(dashed blue line) and the QG model (solid red line). 



We start with the stronger collapse models: the CSL 
and the QG. We will use the d vs. D diagram (Fig. 9) 
to determine how they can be falsified using MERID. 
From the bounds listed in Table I, the decoherence given 
by the collapse model contributes only to conditions (v) 
d < £(t\) and (vii) d < ^3/(At 2 ). Hence, we recalculate 
these two bounds using A = Abb + Acm, where Acm 
is the contribution given by the collapse model; this is 
summarized in Table II. The times ti and t 2 depend 
only on the standard decoherence and are thus chosen 
as described in the caption of Fig. 9. Figure 15 shows 
the d vs. D diagram including either the CSL or the 
QG collapse model. As explained in the caption of the 
figure, the green region is the parameter regime where 
the collapse model can be falsified. This shows that both 
the CSL and the QG collapse model can be tested, for 
instance, at P = 10~ 14 Torr, Ti = 100 K, and \ = 1000, 
for spheres with a diameter ranging from 100 nm to 500 
nm approximately. As expected, this region is larger for 
the QG model than for the CSL model with 7 CSL = 10~ 16 
Hz. The appearance of the green region is mainly due to 
the tighter bound given by condition (vii) d < i/3/(At 2 ) 
which is imposed in order to preserve the visibility of 
the interference pattern. A simulation of the interference 
pattern is also plotted in Fig. 15 for a particular point 
in the diagram, see the figure's caption. The simulation 
is done by numerically solving the master equation at 
the different steps of MERID. Therefore, both the CSL 
(with the conservative value of 7 C sl = 10~ 16 Hz) and the 
QG collapse model can be falsified with the successful 
observation of the interference pattern if MERID were 
implemented at the the green region of the parameter 
regime. 

These results are in strong contrast with the ones ob- 
tained for the much weaker exotic decoherence given by 
the DP and K collapse models. Indeed, the green region 
is zero for both cases even at much higher vacuum con- 
ditions and low bulk temperatures. However, we have 
shown in Fig. 14 that the localization parameter of these 
collapse models is larger than the one given by blackbody 
radiation at Ti = 4.5 K, and thus their additional exotic 
decoherence should be observable. However, the prob- 
lem in this case is that the time required to observe this 
decoherence using MERID is longer than the coherence 
time allowed by the scattering of air molecules 1 /j a i r (see 
Fig. 6). More specifically, the reduction in the visibility 
for d -C 2ocm and t 2 <C 1/ 7air is given by 



V(f a ) = exp [-(Abb + A C M)d%/3] 



(50) 



Assuming Abb "C Acm, the visibility can be approxi- 
mated to V{t 2 ) ~ exp [— AcMd 2 *2/3] , and thus it will be 
reduced at times 



3 / 1 

A C m a" 7 a h- 



(51) 



where we added the second inequality to emphasize that 
t 2 has to be smaller than the coherence time allowed by 
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scattering of air molecules. By inspection of Figs. 6, 
12, and 13, one realizes that these two conditions are 
extremely challenging to be fulfilled for the DP and K 
model. Therefore, we conclude that the DP and the 
K model cannot be tested using the present form of 
MERID. However, the fact that the localization param- 
eter of these collapse models is larger than the one given 
by blackbody radiation at cryogenic bulk temperatures, 
hints at the fact that an improved version of MERID, in 
which, for instance, the free dynamics are accelerated by 
using repulsive potentials, could meet this challenge. 



V. OPTOMECHANICAL DOUBLE SLIT 

The last section of this article is devoted to the analy- 
sis of the restrictions that an optomechanical implemen- 
tation of the squared position measurement imposes on 
MERID. The implementation of MERID using cavity op- 
tomechanics with levitating dielectric spheres has been 
recently proposed in [32]. Here, we provide a thorough 
derivation and we better remark the conditions that need 
to be fulfilled. For further literature on cavity optome- 
chanics with levitating dielectrics, we refer the reader to 
[16, 29, 31, 32]. 

In the following we focus on step 4 of the protocol, 
where the sphere is assumed to enter a small cavity, 
ideally aligned such that the mean position along the 
cavity axis (x) of the sphere is at the node of a cavity 
mode. In this configuration, the optomechanical coupling 
is quadratic with x. This implies that the output light 
of the cavity contains information about x 2 , and there- 
fore, this can be measured by homodyning the light. The 
optomechanical Hamiltonian reads [31] 



H(t) = — + hga)ax 2 + \hE(t)(a - a f ) 
2m 



(52) 



The first term describes the kinetic energy of the sphere 
along the cavity axis (note that that there is no har- 
monic potential since the particle does not need to be 
trapped during the short interaction required to measure 
x 2 ). The third term describes a time dependent driving 
at frequency lul which equals to the cavity resonant fre- 
quency oj c , which is used to parametrize the short light 
pulse. Finally, the second one is the important term de- 
scribing the optomechanical coupling when the sphere is 
placed at the node of the cavity mode. We have defined 
the creation (annihilation) operators of the cavity modes 
(a), the dimensionless position operator x = x/a, and 
the optomechanical coupling rate given by g = go<J 2 /xq, 
where 



fjo 



V 



(53) 



in the case of a nanosphere [16, 29, 31, 32]. Here, e c = 
3Rc [(e r — l)/(er + 2)] depends on the relative dielectric 
constant e r , k c — lu c /c, and V c = TrW 2 L/4 is the cavity 



volume, where W is the waist of the cavity mode and L 
the length of the cavity. As discussed in [32] , note that 
g enhances go by a potentially very large factor a 2 jx\ 
depending on the size of the wave packet. The interaction 
time is assumed to be very small so that the interaction 
is in the regime of pulsed optomechanics [57] . We do not 
take into account the optimization of the pulse shape [57] 
and simply consider a time dependent driving frequency 
given by E(t) = y^27cn^h£(i), where is a flat-top 
function of length T and amplitude ~ 1/vT such that 
Jo £ 2 {t)dt = 1, k is the decay rate of the cavity, and 
n p h is the total number of photons that the light pulse 
carries. The decay rate of the cavity has a contribution 
given by the finesse IF of the empty cavity and by light 
scattering, and it reads [29, 32] 



2tt ce 2 V 2 kl 



2TL 16nV c 



(54) 



A. Measurement operator and strength 

Let us show here that after a short interaction with 
the light pulse, the measurement of the phase quadra- 
ture of the output light realizes a measurement of x 2 
with some given measurement strength x- Pulsed op- 
tomechanics [57] consists in implementing a very short 



interaction time T < 



such that 



(f) T _ (2n + l)uT 
~2m~h ~ 4 



< 1. 



(55) 



This allows us to neglect the kinetic term in Eq. (52), 
which yields 



H(t) ps hga)ax 2 - \%E{t){a ! - a). 



(56) 



In order to obtain the output light quadrature we make 
use of the input-output formalism [58]. The Langcvin 
equation associated to a is given by 



d(i) = -(igx 2 + n)a(t) + E{t) + V2na in (t), (57) 

where ai n is the input cavity noise operator. We fur- 
ther assume that k^> g, such that one can adiabatically 
eliminate the cavity mode by setting a(t) =0. This leads 
to 



o(t) w (E(t) + V2m in (t)) 



1 



• 2 

igx 



(58) 



By using the input output relation a out (t) = \/2na{t) 
ai n (i) and defining the phase quadrature P^it) 
i(aj ut (i) — a out (t))/\/2, one obtains the relation 



PL(t)^Pt(t)+x(t)x 2 , 



(59) 



where x(<) = 2gE(t) / '(k^/k), and we have neglected the 
small term ~ 2gx 2 X\ Ci j n. A balanced homodyne mea- 
surement of the output field performs a quantum mea- 
surement of the time-integrated output quadrature given 
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m a 


Ebb 


8x 


J" 


L 


A 


W c 


2201 Kg/m 3 


2.1 +il0" 10 


2tt x 100 KHz 


0.1 


4.5 K 


28.97 amu 


2.1 + i0.57 


0.1 nm 


1.3 x 10 5 


2 /im 


1064 nm 


1.5 /im 



TABLE III: Experimental parameters used in this article. 



by [32] 



P, L 



4= / Pout(t)dt 



(60) 



An important result is the value of the measurement 
strength, which is given by 



X 



2V2- 



(61) 



An optimization of the pulse shape provides a differ- 
ent pre-factor which slightly increases the measurement 
strength (see [57, 59]). If the measurement of the optical 
phase yields the measurement outcome pl, the measure- 
ment operator describing the collapse of the center-of- 
mass state of the sphere is given by [32, 57, 59] 



A4 = cxp i(f>d s x 



[XL 



XX ) 



(62) 



As a result of the measurement operator of Eq. (62), a 
superposition of two wave packets, separated by a dis- 
tance d = 2 cry 'xl/x an d a width given by approximately 
a d ~ a '/(4\/xlx) = °" 2 /(2dx), is prepared. This is thus 
in full agreement with the treatment of Sec. II C. Fur- 
thermore, the global phase accumulated during the in- 
teraction with the classical part of the held is given by 



g(a\t)a(t))dt 



E\t) 



s 2 (£ 4 ; 



dt 



2grt p h 



(63) 



where the second term can be neglected in the regime 
k > g. 

B. Restrictions to MERID 

The double slit implementation imposes the follow- 
ing restriction on MERID. First, recall that the phase 
(Eq. (63)) needs to be compensated by the phase accu- 
mulated during the time of flight </>tof = t\0j/A, Thus, 
4>tot + 4>ds ~ is fulfilled when the total number of pho- 
tons in the light pulse is 



U)tl k x, 







go c- z 



(64) 



Here we have used again that a 2 ~ x\t\i>j 2 at times t\UJ 
1. Inserting equality Eq. (64) into the definition of the 
measurement strength Eq. (69) leads to 



/, n3/2 go 

V K 



(65) 



Additionally, standard decoherence due to light scat- 
tering during the light-mechanics interaction [29, 32] is 
prevented if the following conditions are met. This deco- 
herence is also of the localization type, with a decoher- 
ence rate given by r sc (t) = A® c E 2 (t)a 2 / k for distances 
smaller than the optical wavelength, where the localiza- 
tion parameter is [29, 32] 



A" 



2;„6 



67T V c 



(66) 



This form of 
J^T sc (t)dt « 
dition on t\ 



decoherence is prevented as long as 
1, which gives rise to the following con- 



;r° 



(67) 



where we have used Eq. (64) and we have defined 
A 



sc A 



Bear in mind that the adiabatic elimination used 
in the derivation is valid as long as k > j = got 2 oj 2 , 
which leads a further constrain, namely t\ <C oj^ 1 \J n/ go- 
Thus, the optomechanical implementation of the step 4 
of MERID yields an additional upper bound on t\ given 
by 



, „ 1 . flT 4g 
h < - mm , /— , -Q- 

w I V 90 r° c 



= t? M . 



(68) 



This is incorporated in Table I as condition (ix). Also, 
by inserting Eq. (68) into Eq. (65) we obtain an upper 
bound for the measurement strength given by 



X <C min < — 
\9o 



1/4 



(69) 



These conditions give a strong limitation for the over- 
all performance of MERID which crucially depend on 
the quality of the optical cavity employed. In [32], 
it was suggested to use the recently developed fiber- 
based Fabry-Perot cavities of length of 2 /mi and finesse 
T w 1.3 x 10 5 [60]. As shown in Fig. 16, for this cavity, 
upper bounds for ii of the order of milliseconds and a cor- 
responding Xmax of several tens are obtained for spheres 
smaller than 100 nm. To see the implications for the 
realization of MERID, we plot in Fig. 17 the d vs. D 
diagram for a pressure of P — 10~ 13 Torr and a bulk 
temperature of Tj = 200 K (which is reasonable consid- 
ering the heating produced by laser absorption [29, 32]). 
For these parameters, t\ = i° M /4 guarantees the ful- 
fillment of the bounds on ti given in Table III. It is 
remarkable that even taking into account the restrictions 
imposed by the optomechanical implementation, spheres 
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FIG. 16: (color online). In the upper panel, ti (see Eq. (68)) 
is plotted as a function of the diameter of the sphere. The 
solid blue line is the upper bound 4go/(ujVg C ) and the dashed 
red line the oj~ 1 (n/go) 1 / 4 , which is the bound given by the 
adiabatic condition. The gray region shows allowed ti values. 
In the lower panel, Xmax (see Eq. (69)) is plotted. We used 
a fiber-based Fabry-Perot optical cavity, see main text and 
Table III. 



with a diameter of tens of nanometers that contain of 
the order of 10 7 atoms can be prepared in superposi- 
tions of the order of their size. Moreover, even the 
CSL model with a localization rate frequency given by 
10 4 

7csL' wn i cn is orders of magnitudes lower than the 
enhancement predicted by Adlcr [20, 48], can be falsi- 
fied. The result shown in Fig. 17 is very similar to the 
one given in [32]. Note however that here we have used 
a pressure three orders of magnitude larger due to the 
saturation effect omitted in [31]. This renders the imple- 
mentation of MERID using cavity optomechanics with 
levitating spheres less challenging. Finally, we remark 
the recent proposal given in [61] to test the CSL model 
using an all-optical time-domain Talbot-Lau interferom- 
eter for clusters with masses exceeding 10 6 amu. 




20 30 40 50 60 



D(nm) 

FIG. 17: (color online). The d vs. D diagram (see Fig. 9) 
is plotted taking into account the additional limitations im- 
posed by an optomechanical implementation of the squared 
position measurement. We used the experimental parameters 
of Table III, a pressure of P = 10 -13 Torr, an internal bulk 
temperature of Ti = 200 K, ti = t? M /4, and t 2 = 10' 2 /j air . 



VI. CONCLUSIONS 

In summary, we have shown that by combining tech- 
niques and insights from quantum-mechanical resonators 
and matter-wave interferometry, one can prepare large 
spatial quantum superpositions of massive objects com- 
parable to their size. The protocol consists of cooling a 
mechanical resonator to its ground state, switching off 
the harmonic potential to let the wave function coher- 
ently expand, preparing a spatial quantum superposi- 
tion by performing a measurement of the squared po- 
sition observable, and observing interference by measur- 
ing the position after further free evolution. We have 
focused on solid spheres with diameters ranging from 
tens of nanometers to few micrometers. We have taken 
into account unavoidable sources of decoherence such as 
the scattering of environmental massive particles and the 
emission, absorption, and scattering of blackbody radi- 
ation. Both sources provide coherence times of the or- 
der of milliseconds within reasonable values for pressures 
and temperatures. At low pressures, decoherence due to 
blackbody radiation is dominant when the bulk temper- 
ature is larger than the cryogenic environmental temper- 
ature. Additional limitations are given by the slow free 
dynamics involved for these massive objects. For larger 
masses, the wave function takes longer to coherently ex- 
pand and to build a visible interference pattern after a 
superposition has been prepared. 

We have also argued that this protocol can be applied 
to test some of the most paradigmatic collapse models. 
In particular, we have analyzed the continuous sponta- 
neous localization model (CSL), a model based on quan- 
tum gravity (QG), the Diosi-Penrose model (DP), and 
the Karolyhazy model (K). The CSL and the DP are 
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much stronger than the DP and K model and can be fal- 
sified using reasonable experimental parameters. In par- 
ticular, the famous continuous spontaneous localization 
model can be tested using the original and conservative 
choice of parameters. However, the DP and K model 
are much more challenging to falsify despite the fact that 
they predict a decoherence which is stronger than the one 
provided by standard decoherence at low bulk tempera- 
tures. Nevertheless, these models are strongly limited by 
the fact that the free dynamics of the large masses re- 
quired is too slow. We remark that for the Diosi-Penrose 
model we did not consider the strong enhancement pro- 
vided by taking into account the mass density at the mi- 
croscopic level. The later, besides being a controversial 
choice, turns the model into a one parameter model given 
by the mass resolution parameter r$. Note that if this 
parameter is taken into account, the protocol proposed 
here provides unprecedented lower bounds to its value. 

We have also addressed the optomechanical implemen- 
tation of the protocol presented, namely MERID. We 
focused on the squared position measurement required 
to perform the double slit, and we have considered cav- 
ity optomechanics with optically levitating nanospheres. 
We have shown that the overall performance of the pro- 
tocol is limited by this implementation, since both the 
global phase added during the interaction and light scat- 
tering set upper bounds on the expansion time and the 
measurement strength. Nevertheless, with recently de- 
veloped fiber-based cavities and for spheres of the or- 



der of tens of nanometers, superpositions of the order 
of their size could be prepared. This provides unprece- 
dented bounds to the continuous spontaneous localiza- 
tion model that can be used to falsify the enhancement 
of the localization rate predicted by Adler. 

There are various directions to further pursue the work 
presented here. First, the study of the implementa- 
tion of MERID using cavity optomechanics with sus- 
pended disks [62]. In this setup, the mechanical fre- 
quency can also be varied since the tight harmonic po- 
tential is achieved by optical trapping, but the scatter- 
ing of light is strongly reduced when the laser waist is 
smaller than the disk. This comes at the price of induc- 
ing decoherence due to the coupling with internal elastic 
modes. Second, the possibility of using repulsive poten- 
tials to exponentially increase the time scales of the free 
dynamics, and thus, to efficiently use the long coherence 
times given by scattering of air molecules and blackbody 
radiation. Naturally, this has to be done without in- 
corporating additional sources of decoherence. In any 
case, we believe that the synergy between the fields of 
quantum-mechanical resonators and matter-wave inter- 
fcromctry will allow to explore in the near future the 
limits of quantum mechanics at unprecedented scales, an 
exciting possibility indeed. 

I am grateful to J. I. Cirac, A. C. Pflanzer, F. Blaser, 
R. Kaltenbaek, N. Kiesel, and M. Aspehneyer for useful 
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